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Abstract. Real linear operators emerge in a range of mathematical 
physics applications. In this paper spectral questions of compact real 
linear operators are addressed. A Lomonosov-type invariant subspace 
theorem for antilinear compact operators is proved. Properties of the 
characteristic polynomial of a finite rank real linear operator are inves- 
tigated. A related numerical function, defined as a normalization of 
the characteristic polynomial, is studied. An extension to trace-class 
operators is discussed. 



1. Introduction 

Real linear operators arise in various applications in mathematical physics, 
such as in Calderon's inverse problem in the plane [2j and in studies related 
with planar elasticity [TTj. An operator i? on a separable complex Hilbert 
space % is said to be real linear if R{x + y) = Rx + Ry and R{rx) = rRx 
for all x,y € "H and r G R. This paper is concerned with compact real linear 
operators and their spectral properties. Compactness means that images of 
bounded sets have compact closures. An invariant subspace result for com- 
pact antilinear operators is proved. To study the spectrum of finite rank real 
linear operators, properties of the characteristic polynomial are investigated 
and a related numerical function is introduced. 

Spectral studies of operators are naturally linked with invariant subspaces. 
For compact complex linear operators, the existence of a nontrivial invari- 
ant subspace was established by von Neumann in an unpublished form. The 
subsequent generalizations in [H [3l O [151 US] ^ead to Lomonosov's theorem. 
We prove an antilinear analogue stating that given a nonzero compact anti- 
linear operator R there exists a nontrivial complex linear subspace invariant 
under every real linear operator in the norm closure of complex polynomials 
in R. Taking into account that a compact real linear operator does not nec- 
essarily have a nontrivial invariant subspace, it is not clear for which class 
of operators the result could be extended further. 

The spectrum of a finite rank real linear operator is given by the zero set 
of a characteristic polynomial p{\, A) with the inclusion of the origin when 
Ti is infinite dimensional. The spectral problem is challenging already in the 
finite dimensional case studied in [UlTJllIl]- The characteristic polynomial 
is a so-called real analytic polynomial so that its zero set, if nonempty, 
is a plane algebraic curve. However, it is nontrivial to show under what 
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assumptions the spectrum is nonempty. For compact real linear operators, 
so far, it is known that the intersection of the spectrum and a straight line 
through the origin consists of isolated points that accumulate at most to the 
origin; see [13] . 

In this paper, it is shown that the characteristic polynomial p can be cast 
in the form 

p{\,\) =v*Hv 

where H = [hij]fj^Q G £;(«,+!) x(n+i) -g ^ Hermitian matrix and v = [X^]'j=o G 
C""*"^. This factorization allows proving a number of properties of p. First, 
p can be represented as a weighted sum of squares where the di 

are real and pi complex analytic polynomials. Moreover, the normalization 

is shown to be connected to the field of values of the coefficient matrix H. 
More precisely, its values are convex combinations of the eigenvalues of H 
and hence a subset of the field of values of H. When transformed to infinite 
dimensions, a similar numerical function for trace-class real linear operators 
is studied by utilizing a generalization of the determinant. 

The contents are as follows. In Section 2 some preliminary definitions are 
given and Shatten classes are discussed. In Section 3 an antilinear version of 
Lomonosov's theorem is proved. In Section 4 the characteristic polynomial 
of finite rank real linear operators is studied. The so-called numerical func- 
tion is introduced and its properties are investigated. A related numerical 
function for trace-class operators is considered. 

2. Real linear compactness 

2.1. Preliminaries. Let Ti he a complex separable Hilbert space. A real 
linear operator on 'H is complex linear if Ri = iR, or antilinear if Ri = 
—iR. With respect to these extremes, a real linear operator R on H can be 
represented uniquely as 

(2.1) R = C + A 

where C = ^{R — iRi) is the complex linear and A = ^{R + iRi) is the 
antilinear part of R. Further, an antilinear operator A can be factored as 
A = Bk, where B = An is complex linear and k is a unitary conjugation, 
i.e., an antilinear operator k that is involutory, = /, and unitary (bijective 
isometry). The complexification of R with respect to k is then a complex 
linear operator on 'H © 'H defined by 

C Ak 



(2.2) R^ 



A real linear operator R is bounded if its operator norm ||i?|| = sup{||i?x|| : 
||x|| = 1} is finite. 

Definition 2.1. The spectrum cr{R) consists of those A G C for which R— X 
does not have a bounded inverse. 

As usual, A is an eigenvalue of R if there exists x / such that Rx = Ax. 
Compactness is defined using (|2.ip . 
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Definition 2.2. A real linear operator R = C + A on T-L is compact if its 
complex linear part C and antilinear part A are compact. 

A real linear operator R is of finite rank n if the so-called left-rank 
n = dim{R'H^)^ is finite. Or equivalently, the dimension n of the smallest 
complex linear subspace of Ti containing the range of R is finite. Obviously, 
a real linear operator is compact if and only if it can be approximated in 
operator norm with finite rank operators, or yet equivalently, if its com- 
plex linear and antilinear parts can be approximated by finite rank complex 
linear and antilinear operators, respectively. 

Fairly recent examples of compact real linear operators are provided by 
the so-called boundary Friedrichs operator and the Friedrichs operator of 
a planar domain. More elaborate applications to inverse problems can be 
found in [2\. 

Example 1. Consider the Hardy space on the unit circle T 

= {/(e^^) = f{k)e'''^ G : />) = for all A: < 0} 

with the inner product {f,g) = ^ f^"" /(e^^)g(e^^) dO. The so-called bound- 
ary Friedrichs operator Ba with symbol a = J2kez O'kS.^^^ G L°°(T) is defined 
on as 

Baf = Paf, 

where P is the orthogonal projection from onto H^. For the orthonormal 
basis {e^'^^I^Q we have 

Thus Ba is an antilinear self-adjoint operator on and the infinite matrix 
with respect to {e^'^^j^Q is a Hankel matrix. Furthermore, it is known 
that Ba is compact when a G C(T). For further reference on more general 
Friedrichs operators, see [171 H] • 

Example 2. Let $7 be a planar domain. The Friedrichs operator with symbol 
a G L°°(r2) is defined on the space of analytic square- intergable functions 
AL'^{n) as 

Fa : AL\Q) ^ AL\Q), Faf = PaJ, 

where P is the orthogonal projection from onto AL'^{Q), called the 

Bergman projection. It is known that if boundary of O is smooth enough 
and if a G C(r2), then Fa is compact. Here O denotes the closure of fl. 
Again, for a more accurate description see pTj 118]. 

Recall that a compact complex linear operator is in the Shatten p-class 
if its singular numbers belong to P. The singular numbers of a compact 
antilinear operator A are defined identically as Sj{A*A), and if they belong 
to P, A is in the Schatten p-class. Then we set the following for compact 
real linear operators. 

Definition 2.3. A compact real linear operator R = C + A on T-L is in the 

Schatten p-class, denoted by Bp{'H), for 1 < p < oo, if both its complex 
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linear part C and its antilinear part A are in Schatten p-class. The Schatten 
p-norm of R is defined as tlie sum 

(2.3) \\R\\p = \\C\\p+\\A\\p 

witli \\Cfp = Y.T=iSj{Cy and \\A\\l = Y.T=iSj{AY. Following standard 
nomenclature, operators in BiiT-L) are called trace-class (or nuclear), and 
those in B2{Ti) Hilbert-Schmidt. 

This is just one way of defining the Schatten norm of a real linear operator. 
It is not clear how different approximations numbers in general are related 
to this definition of singular numbers for real linear operators; see [TTj for 
real linear operators of finite rank. This choice of Schatten norm for a real 
linear operator is equivalent to its complexification in the following sense. 

Proposition 2.4. Let R = C + A be a compact real linear operator on 
Ti with complex linear part C and antilinear part A, and let R^ he the 
complexification with respect to some unitary conjugation k. Then R^ G 
Bp{7i © v.) if and only if R & Bp{H). In addition, it holds 



||(C + ^)^||^>||C||^+P||^ 



and 



l/~tC IIP 

Hp 



1^ Hp 



Proof. Assume R^ & Bp{'H ® K) ■ We have 



I-"- Hp 



sup < ^ |(ii^/ij,/ij)w®w| : {/ij} orthonormal > . 



Choosing {hj} to be {(ej, 0), (0, /j)} where {cj} and {/,} are orthonormal 
sets of we get 

\{R%.h,)u^-H\' = \iCe,.e,)\^ + Y UKf„f,)\P. 
j j j 

Since the orthonormal sets {ej} and {/,} are arbitrary, we see that C and 
Ak, hence and also A, are in Bp{l-L). 

Assume R = C + A & Bp(n) so that by definition C,A £ BpCH). It 
follows for the complexifications C'^,A'^ G Bp{'H (BTi). Namely, as we have 
the representation C = J2j ^j)fj some orthonormal bases {cj} 

and {/}■}, then kCk = Sj{C){-, Kej)Kfj. Then for all {x,y) £ V. (Bl-L it 
holds 



'c 




X 


kCk 




y. 



X 






) 




y 



















Kfj^ 



where the sets {(cj, 0), (0, nej)} and {{fj, 0), (0, nfj)} form orthonormal bases 
oi H (B Ti- This means that the singular values of are those of C but 
each repeated twice and therefore ||C"^||p = 2||C||p. 

In a similar fashion, the representations Ak = Sj{A){-, Kej)fj and 
kA = J2j ■5i(^)('i ^j)i'^fj lead to 



Ak 



kA 







+ Sj{A) 
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Again, it is seen that \\A^\\f, = 2\\A\\f,. It follows G BpiV. U). □ 

Although the trace-class of real linear operators has been defined above, 
defining an actual trace for a real linear operator is seemingly more problem- 
atic. Namely, for an antilinear operator A, the sum Y^ji^^j^^j) depends 
on the orthonormal basis {ej}. Fixing an orthonormal basis of Ti allows 
defining a basis-dependent trace of a real linear operator as the trace of its 
complexification with respect to this basis. 

3. Invariant subspaces 

As is well known, eigenspaces corresponding to specific eigenvalues deter- 
mine invariant subspaces for a compact complex linear operator. Recall that 
a set C is invariant under an operator L when LA4 C M. Granted, 
a specific eigenvalue determines a real linear invariant subspace for a real 
linear operator. However, as there are real linear operators defined more 
naturally on complex linear Hilbert spaces, studying complex linear invari- 
ant subspaces is, for the same reason, more natural. This is also done to 
preserve the structure of the underlying complex Hilbert space. 

In general, for a real linear operator, having a complex linear invariant 
subspace is quite a restricting property. 

Proposition 3.1. Let M he a complex linear invariant subspace under a 
real linear operator R = C + A onH, where C and A are the complex linear 
and antilinear parts of R. Then Ai is invariant under both C and A. 

Proof. Assume x € A4. Then also Cx = ^{R — iRi)x € Ad and Ax = 
\{R + iRi)x ^ M. □ 

For this reason, we are concerned with compact antilinear operators and 
show that then there exists, much like in the complex linear nontrivial 
closed complex linear invariant subspace. The investigation of the invari- 
ant subspace problem is an longstanding project. In 1930, von Neumann 
proved that every compact operator has a nontrivial closed invariant sub- 
space. Aronszajn and Smith extended this to compact operators on Banach 
spaces [T]. Using non-standard analysis, Bernstein and Robinson proved 
the version for polynomially compact operators on a Hilbert space [3J, and 
Halmos proved the same using classical analysis [8J . Lomonosov proved that 
there exists a nontrivial closed subspace of a Banach space that is invariant 
under every operator commuting with a given nonzero compact operator 
|15j . Hilden simplified the proof of this [16] . For recent advances and a 
modern approach on the invariant subspace problem in general, see [3]. 

Let p{z) = Yl^=o'^j-^'^ be a polynomial with aj G C. A complex poly- 
nomial in a real linear operator R is defined to be the real linear operator 

Definition 3.2. For a compact antilinear operator A on H, define {A) 
to be the closure in the operator norm of polynomials in A with complex 
coefficients. 

Note that {A) is almost a subalgebra in the sense that if A G C and 
T,S £ (A), then AT, T + S, TS G {A). Obviously the main obstruction 
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is T{\S) 7^ ATS". In general, the elements of (A) are genuinely real linear 
operators since even powers of A are complex linear while odd powers are 
antilinear. Also, (A) contains noncompact elements, e.g. all polynomials 
p{A) = CLj^'' foi' which ao 7^ 0. Clearly 7W is a closed complex linear 
invariant subspace under A if and only if it is such under every T G {A) . 

Like in the complex linear case, invariant complex linear subspaces for 
an antilinear operator A may be generated by forming the closure of the 
complex linear span of {A^y}j?^Q for any nonzero y (z "H. Indeed, denoting 
by span£ the closed complex linear span, 

(3.1) span{A^y}~o = {A)y = {Ty : T G {A)}, 

c 

is a closed complex linear subspace of Ti. Either {A)y is a closed nontrivial 
invariant subspace under A, or {A)y = % whence y is called a cyclic vector 
of A. 

The arguments of Lomonosov and Hilden |16] largely apply in proving the 
following. Of course, antilinearity requires some changes. Indeed, the fact 
that is a complex linear subspace of % is crucial. 

Theorem 3.3. Let A he a nonzero compact antilinear operator on %. Then 
there exists a nontrivial closed complex linear subspace ofH that is invariant 
under every operator T £ (A) . 

Proof. Assume the contrary, that there is no such subspace. If A G C is an 
eigenvalue of A with the corresponding eigenvector x & Ti, then spandx} 
provides a contradiction. Thus it must be that a{A) = {0}. Without loss 
of generality, we may assume ||^|| = 1. 

Let xq € H he such that H^xqH > 1. Then we have 

(3.2) ||xo|| = Pllllxoll > ll^a^oll > 1 

from which it follows that ^ i?i(xo), the closed unit ball with centre at xq. 
Moreover, we have ^(-61(2:0)) C Bi(Axo), the closed unit ball with centre 
at Axq, because for all x G Bi{xo) it holds \\Ax — Axo\\ < \\A\\ \\x — xo\\ < 1. 

Thus, ^ A(Bi{xo)) since \\Axo\\ > 1. 
For all nonzero y £ H, 

Vy = {p{A)y : p{A) is a complex polynomial in A} 

is a complex linear subspace invariant under every polynomial p{A) G (A). 
Thus its closure is a complex linear closed subspace invariant under every 
T G (A). As Vy {0} it is dense in V. due to the counter-assumption. It 
follows that for every nonzero y £ T-L there exists a polynomial p{A) G {A) 
such that ||p(A)y — xo|| < 1. For each polynomial p{A) G (A), define the 
open set f^p(yi) = {y € H : \\p{A)y — xq\\ < 1}. The collection {C/p(^) : 
p{A) is a polynomial in A} is an open cover of ?{\{0} and thus of A{Bi{xo)) 
since ^ A{Bi{xo)). By the compactness of A, A{Bi{xo)) is compact, and 
there is a finite subset Vq C (A) of polynomials such that {C/p(yi) : Pi^) ^ "^o} 

covers A{Bi{xo)). Consequently, for every y G A{Bi{xq)), there exists a 
p{A) G Vo such that \\p{A)y - xo\\ < 1. 

Since Axq G ABi{xo), there exists pi{A) G Vq such that |[pi(A)j4xo — 
xoll < 1. Thus Api{A)Axo G ABi{xo), and there exists P2{A) G Vq such 
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that \\p2{A)Api{A)Axo — a;o|| < 1- By induction, we can define a sequence 
{pn{A)} in Vq such that 

(3.3) \\pn{A)A - ■ ■ p2{A)Api{A)Axo - xoW < 1. 

For a polynomial p{A) = Ylj=o '^i^"') the polynomial with conjugated co- 
efficients p{A) = Ylf=o unique polynomial that satisfies Ap{A) = 
p{A)A. Suppose there are two polynomials pi{A) and P2iA) for which 
Ap{A) = pi{A)A = P2{A)A. Then q{A) = 0, where q{A) = pi{A)A - 
P2{A)A = '^fJi'bjA^. But the coefficients of q{A) must be all zero, for if 
the leading coefficient 6fe 7^ 0, 1 < /c < + 1, then span^{A^ z}^~q with an 
arbitrary 2; 7^ is a complex nontrivial closed subspace under every T € {A) . 
This is in contradiction with the counter-assumption. 

Let M = max{||p(^)||, ||^(^)|| : p{A) £Vo}. Then 

(3.4) \\pr,{A)A---p2{A)ApiiA)A\\ = ||p„(A) • • •pi(A)^"|| < ||(MA)"||, 
for all n = 1, 2, . . ., where Pk{A) is either PkiA) or pj^{A). 

Lemma 3.4. It holds that lim„^oo IKMA)*^!] = 0. 

Proof. The complex linear compact operator (MA)'^ has no eigenvalues, and 
the orgin is the only spectral value. For if it had an eigenvalue A with the 
eigenvector x, then span^-jx, Ax} would be a nontrivial complex linear in- 
variant subspace, contrary to the counter-assumption. Thus, by the spec- 
tral radius formula, it holds that lim„_!.oo IK-M^)^"!!"^/" = which implies in 
particular that lim^^oo ||(MA)2"|| = 0. Then also lim„^oo ||(MA)2"+i|| = 
which proves the claim. □ 

Using (fXi]) with Lemma E31 we get lim„^oo ||Pn(^)^- • -^1(^)^11 = 
which combined to (j3.3|) leads to ||a;o|| ^ 1) ^ contradiction to (|3.2|) . thus 
proving the antilinear Lomonosov's theorem. □ 

It is not clear how Theorem 13.31 could be strengthened. In the proof 
above it is crucial that {A)y yields a complex linear subspace. In general, 
for a compact real linear operator R a similarly defined set {R)y cannot 
be expected to be invariant under R. In particular, a compact real linear 
operator may not have any nontrivial complex linear invariant subspaces. 
An example of this is the operator 



"1 


1' 




"0 


0" 





1 


+ 


1 






on C^, where r denotes the componentwise complex conjugation. 

Finally, let us comment on the connection between eigenvalues and finite 
dimensional invariant subspaces. In the complex linear case, the existence of 
a finite dimensional invariant subspace entails the existence of eigenvalues. 
This reasoning is not valid for real linear operators. An obstruction is given 
by the following proposition, where the minimum modulus (or injectivity 
modulus) is defined as j{R) = inf{|[i?x|| : ||x|| = 1}. 

Proposition 3.5 ([14J). Let a real linear operator R = C + A be split as 
R = R + T , where T = ^{A — A*) is the skew- adjoint part of the antilinear 
part of R and R = C-'r^{A+A*). < j{T), then R has no eigenvalues. 
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Proof. By the antilinear skew-adjointness of T, Tx is always orthogonal to x 
for any x gH. Assume that A is an eigenvalue of R with eigenvector x 
of unit length. Then we have by the Pythagorean theorem 

|A|2 + \\Txf = \\Xx - Txf = \\Rxf. 

Taking infimum on the left-hand side and supremum on the right-hand side 
over unit vectors, we have 

|A|2+j(r)2 < j|i?||2 

which contradicts ||i2j| < j(T). □ 

When H is infinite dimensional, the origin is always in the spectrum 
of a compact real linear operator. However, all points in the spectrum 
of a compact real linear operator outside the origin are eigenvalues. In 
this regard, to determine whether there is an associated invariant subspace 
Proposition 13.51 is crucial. Indeed, assume the Hilbert space splits into the 
direct sum T-L = Hj , where each Tij is a finite dimensional invariant 

subspace for the compact real linear operator R. If the assumption \\Rj\\ < 
j{Tj) holds for the restrictions Rj of R to Hj, then R has no eigenvalues. 
In particular, from this it can be seen that cr{R) = {0} does not imply 
lim„_>.oo = for real linear operators in general. 

For this reason, the eigenvalue problem for compact real linear operators 
is challenging in general. Also, it underlines the spectral questions for finite 
rank real linear operators, which we will discuss next. 

4. Characteristic polynomials and numerical functions 

The real linear eigenvalue problem, i.e., characterizing the points A € C 
for which i? — A is singular, is a nontrivial task already in finite dimensions. 
We treat the finite dimensional case first in terms of the characteristic poly- 
nomial of R and a so-called numerical function. Then we discuss a slightly 
different numerical function for real linear trace-class operators which is 
defined through a limit process of finite rank operators. 

Finite rank operators are best dealt with real linear matrix analysis. I.e., 
we assume having a represention of R on as 

(4.1) R = C + Bt 

with C,B ^ C""^". Here r is the componentwise complex conjugation tz = 'z 
for z G C". This representation is attained with respect to some orthonormal 
basis {cj}"^]^ of the complex linear span of RTi. 

Definition 4.1. The characteristic polynomial of a real linear operator R 
on C" is defined as the polynomial 

(4.2) p{\ A) = det(i? - A)^ = det 

The characteristic polynomial is not a complex analytic polynomial but 
a so-called real analytic polynomial. Calling it characteristic is justified by 
the fact that 



C-\ B 
B C -X 



(4.3) 



a{R) = {A G C : p(A, A) = 0}. 
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Polynomials defined in this way have been studied in [6j and [12j , where they 
are called characteristic bivariate polynomials resulting from the separation 
of X = a + iP into its real and imaginary parts. 

The characteristic polynomial can be cast in a canonical form involving 
complex analytic polynomials. 

Theorem 4.2. The characteristic polynomial p{X, A) = Y17j=o ^ij^^'^ of a 
real linear operator R on C"" satisfies the following: 

(i) p is real valued. 

(a) The coefficient matrix H = [/ijjJ^^^Q is Hermitian. 
(Hi) Let U = [uij]lj^Q £ c(n+i)x(n+i) unitary such that UHU* is diag- 
onal. Then p can be cast in the form 



(4.4) p{X,X) = J2d^\p^{X)\^ 

where G R and 



1=0 



PiW = fJ-ij^^ for i = 0,l, 
j=0 



,n. 



The polynomials pi have no common zeros, 
(iv) If H is positive definite, then p can be cast in the form 



p(a,a) = ^Ik(a)P 



i=0 

where each pi is a complex polynomial of degree i. The polynomials pi 
have no common zeros, and consequently, cr{R) = 0. 

(v) IfdetR^ < 0, thenp{X,X) has a zero, and consequently, cr{R) ^ 0. 

(vi) For the adjoint R* , it holds pj?*(A,A) = p(A,A), i.e., the coefficients 
of PR* are complex conjugates of p. In particular, the characteristic 
polynomial of a self- adjoint operator has real coefficients. 

Proof. By properties of the determinant, the characteristic polynomial is of 
degree 2n and can be represented as 



(4.5) 



p{X, A) = ^ hijX^X' = vlHvx, 

i,j=0 



where vx = (1, A, A^, . . . , A") G C"+i and H = [hij]^^^^^ G 
the coefficient matrix. 



(i) Since 



p{X, A) = det 
= det 



C-X B 
B 

/" " 
/ 



C-X B 
B 



/ 
/ 



(n+l)x(n+l) 



IS 



p{X,X), 



p is real valued. 
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(ii) Since p{X, A) is real valued, the two polynomials 

n 

i,j=0 

and 



p(A,A) = h^jX^X' = hjiX^^ 

i,j=0 i,j=0 

are identical for all A € C. Then the coefficients corresponding to 
identical powers A-' A* must also be identical, i.e., hij = hji. 
(iii) As the coefficient matrix is Hermitian, we have H = U*DU with U = 
[uij] unitary and D diagonal with real diagonal entries di. Then 



p{X,X) = vlHvx = (UvxTDUvx = ^(i,|p,(A)| 



i=0 

where Pi{X) = Ylj=o'^ij^'' ■ '^^^ polynomials cannot have common 
zeros since, if there exists /i G C such that = for alH = 0, . . . , n, 
then Uv^ = contradicting the unitarity of U. 

(iv) Let P be the order reversing permutation matrix {xq,xi, . . . ,Xn) i— 
{xn, ■ ■ ■ ,xi,xq). Then p{X,X) = u\P* HPu\, where u\ = Pv\ = 
(A"", • • • , A, 1). The matrix H = P*HP is Hermitian and positive defi- 
nite so that it has a Cholesky decomposition as H = T*T, where T is 
upper triangular with strictly positive diagonal elements. Then 

n 

p{X,X) = ulT*Tux = {TuxTTux = Y \P^W\^^ 

i=0 

where pi is of degree i. The existence of a zero common to all of the 
polynomials pi would contradict the invertibility of T. Then p(X, A) > 
for all A E C, or equivalently (7{R) = 0. 

(v) The leading term of p{X, A) is |Ap". Restricted to the real line /(r) = 
p{r, r) is a continuous function with /(O) = det R'^ < and /(r) — > +oo 
as r — > oo which proves the claim. 

(vi) By the fact that taking the transpose leaves the determinant unchanged 
and the characteristic polynomials are real valued, it holds 



p{X, A) = det 



C-A _B 
B C-X 



det 



'C*-X 

B* C* -X 



= PR*iX,X). 

Denoting the coefficient matrix oi pR* by [vij], it follows that 



PR*{X,X) = Y njX^X'=p{XA) = Y ^^1^'^' 

i,j=0 i,j=0 

for all A G C. Thus rij = hij. Evidently then, for self-adjoint R, the 
coefficients are real. □ 

Theorem 14.21 above connects the characteristic polynomials of real linear 
operators to the study of positive polynomials and sums of squares which is 
classical. Sum of squares decompositions arise in classical moment problems. 
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real algebraic geometry, optimization, and control theory problems. See [9] 
for a thorough exposition. 

The coefficient matrix H cannot be negative definite since the leading 
coefficient is hnn = e.^Hen = 1, where = (0, . . . , 0, 1). Indefiniteness of H 
is a necessary condition for the nonemptiness of (t(-R)\{0} but not sufficient. 
An example of such an operator is 



(4.6) 



A 
A 



with A 



ll+e 




2 




ll-e 


ll+e 


/ 2 


V 2 



and < e < 1. 



The characteristic polynomial is then p(A, A) = |A|^ — 2e|Ap + 1 for which 
the coefficient matrix is diagonal with entries 1, — 2e, 1. 

Recall that the field of values (or synonymously numerical range) of a 
matrix A G C"^" is the closed convex subset of the complex plane 

W{A) = {x*Ax : \\x\\ = 1}. 

For the field of values, see, e.g. [TD]- By being the closed interval between 
the smallest and largest eigenvalue of H, the field of values of a Hermitian 
matrix H is easily computed. 

In terms of the field of values, Theorem 14.21 clarifies that, strictly speak- 
ing, any given Hermitian matrix does not give rise to the characteristic 
polynomial of a finite rank real linear operator. Furthermore, on C the 
real linear operator a + /3t with a, /3 € C has the characteristic polynomial 
p{X, A) = |ap — |/3p — ZSA — aA + AA. Thus, there are Hermitian matrices that 
are not scalable or translatable to the coefficient matrix of a characteristic 
polynomial. 

To connect further the field of values of the coefficient matrix H to the 
characteristic polynomial p of a real linear operator R, the following nor- 
malization is introduced. 

Definition 4.3. Let Rhea real linear operator on with the characteristic 
polynomial p. Then 

p(A,A) _vlHvx 



F : C 



F{X) 



is said to be the numerical function of R, where = (1, A, . . . , A"). 

Based on the properties of the characteristic polynomial, we have the 
following. 

Theorem 4.4. Let F : C be the numerical function of a real linear 

operator R on C", and let the coefficient matrix of the characteristic polyno- 
mial be diagonalized as H = U*DU with D diagonal and U unitary. Then 
the following hold: 

(i) F is continuous and its range is a subset ofW{H). 
(a) F{X) = if and only if p{X, A) = 0. 

(Hi) For every X £ C, F{X) is a convex combination of eigenvalues of H, 

b.(A)P 



(4.7) 



F{X) = Y,d^ 



1=0 
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where the polynomials pi and the real numbers di are as in Theorem 

[23 (Hi). 

(iv) In particular, 

n 

F{0) = detR'^ = ^di\uio\^, 

i=0 

n 

lim F{X) = y^di\uin\'^ = I. 

so that 

[deti?^,l) C F(C) C WiH) = [Xn.in{H),\^,.{H)], 

where Amm(-f^) ^.n-d Amax(^) o-i^^ the smallest and greatest eigenvalue 
ofH. 

Proof. (i) From the form of F it is seen that F is continuous from C to 
R. For any A G C 

(ii) As (X^"=o l-^P'')^^ > for ah A G C, p vanishes exactly when F van- 
ishes. 

(iii) Let p(A, A) = Yl^=o ^ilP'i-Wl'^ representation of Theorem 
Ol (iii). Since U is unitary, there holds ^^^=0 I^P^ = EHo IPiWl"^- 
Thus 

Ik(a)P 



which is a convex combination of eigenvalues {di} of i?. 
(iv) From the definition of F it is seen that F{0) = det R!^ and from Equa- 
tion (gT]) that 

n 

F{0) = Y,di\uio\^. 

i=0 

It holds 

J. _p(A) I A + (lower order terms) ^ 

^"=0 lAp-' lAp*^ + (lower order terms) 

On the other hand •rin''^^}\,ij — )• |umP as |A| — )• oo so that 

n 

lim F(A) = y^di\uin\'^, 

which proves the equality. 

As F is continuous and H is Hermitian, we have the following in- 
clusions: 

[deti?^,l) C F(C) C W(H) = [Xmin{H),Xm.AH)]. □ 
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The numerical function F can be defined for an arbitrary Hermitian ma- 
trix H or even for any matrix A. As exemplified by ()4.6p . the closure of 
F{H) need not equal W{H). It is of independent interest to inquire how 
much of the field of values W{H) remains uncovered by F{C) in each case. 
Numerical evidence suggests that some always remain. The interest stems 
from the fact that calculating the range of F[C) is computationally a one 
dimensional problem whereas the field of values is defined on the unit sphere 
of C^. 

In the subsequent defintion of the numerical function for trace-class op- 
erators, we will need the following rotation property. 

Proposition 4.5. Let F , F^q he the numerical functions and p, p-o the 

characteristic polynomials of finite rank real linear operators C + Bt and 
e-*^C + Bt, where 6* € M. Then it holds for aUr>0 

F{e'^r) = F^0{r) and p{e'^r, e''^r) = p^eir, r). 



Proof. By properties of the determinant, 



P-e(r, r) = det 



ie/2j- 



ie/2j 



-10 



C 
B 



B 

e'^C - 



ie/2j 



i0/2j 



p{e''"r, e 



As v\v\ with A = e r is indepenent of the argument e 
for the numerical function as well. 



ie 



the identity holds 

□ 



Recall that for a complex linear trace-class operator C the determinant 
can be extended in the following way 

dei{I + C)= lim det(/-FC„), 

where \\Cn — C\\i — )• as n — )• cxd. Furthermore, the function / : C ^ C, 
f{z) = det(/ -|- zC) is entire and depends continuously on C in the sense 

|det(/ + Ci) -det(/ + C2)| < ||C7i -C72||iexp(l + ||C7i||i + ||C72||i), 

where Ci and C2 are complex linear trace-class operators. Also, / -|- C is 
boundedly invertible if and only if det(I -|- C) / 0. See [5l Section XI. 9] or 
[3 Chapter II] for details. This allows to define the following. 

Definition 4.6. Let R = C + Ahe a real linear trace-class operator on %. 
Define its characteristic function for A = re*^ 7^ as 



ip{re 



10 \ 



det 



,-1 



{e-'^C + A) 



Proposition 4.7. The characteristic function (p of a real linear trace-class 
operator on H is continuous on C \ {0} and its zeros are the eigenvalues of 
R. Furthermore, for all X ^ 



lim 



PniX,X) 



lA 



2n 



where Rn = Cn + An are finite rank operators with characteristic polynomials 
Pn{X,X) such that — R\\i — > 0. 
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Proof. For fixed 6* € M, we have ip{re^^) = fe{^), where fg{z) = det(I — 
z{e~^^C + A)*^) is an entire function. Thus restricted to the half-line {re*^ : 
r > 0}, If is continuous. Furthermore, by the properties of f{z), f{r^^) = 
with a positive real if and only if r is an eigenvalue of (e~*^C + A)'^. It 
is known [13J that M n a(e-*^C + ^) = M n cT((e-^^C7 + Af). 
On the other hand, for fixed r > 0, we have that 

< \e'^ - e^'^IIICIIi exp(l + \\{e-'^C + Af\\i + IKe-^'^C + ^)^||i). 

Thus is continuous in the angular direction as well, and continuous in 
C\{0}. 

For the limit, we have for all A = re*^ € C \ {0} 



= det(/ - r-'{e-''Cn + A^))^ ^ f{re'') 

as n — 7> oo by properties of the 2n-by-2n determinant defining the charac- 
teristic polynomial. □ 

In Proposition 14.71 use has been made of Proposition 12.41 in the equiva- 
lence of convergence of real linear operators and their complexifications in 
Schatten norms. 
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